Kosterlitz-Thouless and Manning Condensation 
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A comparison between the Kosterlitz-Thouless theory of metal insulator transition in a two 
dimensional plasma and a counterion condensation in a polyelectrolyte solution is made. It is 
demonstrated that, unlike some of the recent suggestions, the counterion condensation and the 
Kosterlitz-Thouless transition are distinct. 



I. INTRODUCTION 

The polyelectrolyte solutions present one of the out- 
standing problems in the field of physical chemistry. Al- 
though the work of Debye and Huckel [ffl shed some light 
on the unusual behavior of simple symmetric electrolytes, 
the properties of the polyelectrolyte solutions still remain 
to large extent ununderstood. In general a polyion can 
be a flexible polymer chain or a rigid molecule of spe- 
cific shape some of whose monomers are ionized. For 
now, we shall confine our discussion to rigid polyelec- 
trolytes whose polyions can be modeled as cylinders or 
spheres, since as the reader will see, even these signifi- 
cantly simpler systems already present more than enough 
complications. In particular even such basic question, 
as what is the form of electrostatic interaction between 
two polyions, still remains controversial. A prototypical 
molecule whose shape can be modeled as a cylinder is a 
DNA. It should be noted, however, that even in this case 
the approximation of replacing a long molecule by an in- 
finite cylinder will only be valid if the persistence lenght 
is sufficiently greater than the Debye screening lenght in- 
side the solution. As an example of spherical polyions, 
one can consider various colloidal suspensions, for exam- 
ple the ones made of latex or polystyrene particles. 

In this note we shall confine our attention to the case 
of rigid cylindrical polyions. In particular, our discussion 
is motivated by the recent observation of Kholodenko 
and Beyerlein m that the commonly observed counterion 
condensation during which a finite fraction of disassoci- 
ated counterions recondenses onto polyions, is a "special 
case of the Kosterlitz-Thouless || (KT)-like phase transi- 
tion" 0] . Here we shall demonstrate that the KT transi- 
tion and the Manning's counterion condensation Q] have 
nothing in common. The erroneous conclusion reached 
by the above authors can be attributed to the confusion 
over the system of units used in the original KT theory. 

To begin our discussion let us first define the Primi- 
tive Model of Polyelectrolyte (PMP) [|. It consists of 
N cylindrical polyions of a cross-sectional diameter a p 
and lenght L, each carrying P ionized groups of charge 
—q spaced uniformly, with a separation b, along the axis 
of a cylinder. A total of PN counterions are present 
to preserve the overall charge neutrality of the system. 
The counterions are assumed to be spherical in shape 
with a diameter a c , each earring a charge q. The whole 



system is confined inside a region of volume V. It will 
prove convenient to define the distance of closest ap- 
proach between the centers of a polyion and a counterion, 
a = (a p + a c )/2. The solvent is modeled as a uniform 
medium of a dielectric constant D. The bare interaction 
between a sufficiently long polyion and a counterion can 
then be approximated by 



0(r) 



-2q(a/D)ln(r/a), 



r < a, 
r > a, 



(1) 



where a = —Pq/L = —q/b. If we make the replace- 
ment a — > —q (or b — > 1), the logarithmic potential 
which appears in the above expression is exactly the 
same as for two oppositely charged hard disks in two di- 
mensions (assuming we keep the same Poisson equation, 
V 2 </> = —A-Kp q /D, as in three dimensions). Furthermore, 
it is well known that when the temperature of a two di- 
mensional plasma is lowered below a certain value, all the 
ions associate into dipolar pairs, changing the properties 
of plasma from being a conductor to an insulator ||. 
This is a continuous thermodynamic transition, which 
happens to be of infinite order. Is it possible that the 
counterion condensation observed in polyelectrolytes is 
a realization of this transition? In order to answer this 
question we shall present two simple mean- field theories, 
one for a two dimensional Coulomb gas and another for 
a polyelectrolyte solution in three dimensions. 



II. 2D COULOMB GAS 

Our system will consist of hard disks of diameter a, 
half of which carry the charge +q, while the other half 
carry the charge — qF. As is usual for the restricted prim- 
itive model (RPM), the solvent is modeled as a uniform 
medium of a dielectric constant D. Due to strong elec- 
trostatic interactions we expect that some of the ions 
will associate forming dipolar pairs |6). The total den- 
sity of hard spheres is p = p+ + p_ + 2p2, where 
p + = p- = p\/2 is the density of free unassociated ions, 
while the P2 is the density of dipolar pairs. In order to 
explore the thermodynamic properties of this plasma, we 
need the free energy for the system. We shall construct 
this free energy out of two parts: the entropic free energy 
associated with the momentum degrees of freedom in the 
partition function, and the electrostatic free energy due 
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to the interactions between the ions. It is a simple matter 
to write the expression for the entropic part of the free 
energy, since it corresponds to that of an ideal gas. Defin- 
ing the Helmholtz free energy density as / = — F/V we 
have 



E 



Ps [1 - ln( Ps Af l/C.)], 



(2) 



where d is the dimensionality, (3 — l/fc^T, A s is the de 
Broglie thermal wavelength, £ s is the internal partition 
function for species s = +,—,2, and | s | is how many 
particles constitute a specie, i.e. | + | = | — | = 1, | 2 | = 
2. The electrostatic free energy is easily calculated in the 
spirit of the original Debye-Huckel theory M. We find 
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ln[naKi(Ka) 



(3) 



where Ki is the modified Bessel function of first order, 
(ko) 2 = Anpl/T*, and the reduced density and temper- 
ature are respectively p* = pa 2 and T* = k B T/q 2 . The 
densities of the free and the associated ions are related 
through the law of mass action p,+ = /i_ = where 
the chemical potentials are determined from the free en- 
ergy fi s = -df/dp s . We find 



P2 = \K{T)p\e^ 



where the excess chemical potential is pf 



(4) 

-dr i /d Pl 

and the equilibrium constant K{T) — (^(TVC+C-- F° r 
the purpose of the present exposition it is sufficient to 
know that the equilibrium constant remains finite in the 
limit pi — ► 0, for further discussion we refer the in- 
terested reader to ||. In the limit of small densities 
the excess chemical potential can be expanded in pow- 
ers of Pl ,p ex = -[j E + ln( K a/2)]/T* + 0{ Pl ). Sub- 
stituting this expression into Eq. 3, we observe that 



P2 



pfe 



■Hp"i)/t* 



Pi 



*(2-l/T*) 




In particu- 



lar, we see that for a fixed density p, as the temperature 



is reduced to T* — > T^ T = 1/2 the density of free ions, 
pi — ► 0. This is exactly the Kosterlitz-Thouless metal 
insulator transition. Since the Debye screening length is 
inversely proportional to y/pi, we see that it will diverge 
as T ► [L > zr r r • 

Let us now take a look at what happens in the case of 
cylindrical polyelectrolytes || . 



III. POLYELECTROLYTE SOLUTION 

We shall work in the context of the PMP defined above. 
Just as in the case of simple electrolyte addressed in the 
previous section, we expect that the strong electrostatic 
interaction between the polyions and the counterions will 
result in formation of clusters composed of one polyion 
and 1 < n < Z condensed counterions. We are then lead 
to two conservation laws: 
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n=0 



pp = p+ + E npr ' 

n=Q 



(5) 



(6) 



where p is the total density of the polyions, pa is the den- 
sity of free polyions, p n is the density of the clusters of 
size n, and p + is the density of free, unassociated, counte- 
rions. In order to explore the thermodynamic properties 
of the PMP we require the free energy. As in the case 
of simple electrolyte, this can be constructed as a sum of 
the entropic and the electrostatic contributions. The en- 
tropic part of the free energy is given by the Eq. || with 
s = {+,0 < n < Z}. The electrostatic contribution, 
however, is significantly more complex than for simple 
electrolyte. It consists of the free energies due to the 
polyion-counterion, polyion-polyion, and the counterion- 
counterion interactions. In general it is quite difficult to 
take a full account of all these effects. The condensation 
phenomena that we are interested in studing, however, 
already occurs in the limit of vanishingly small densi- 
ties, p — > 0. In this limit, it is possible to show that 
the polyion-polyion and the counterion-counterion inter- 
actions are small and the main contribution to the elec- 
trostatic free energy comes from the polyion-counterion 
interaction This contribution can be calculated in 
the spirit of the Debye-Huckel theory and we find || 
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where (na) 2 = Anp* + /T* and the reduced density and 
temperature are now p* = pa 3 and T* = kBTDa/q 2 re- 
spectively. The distribution of cluster sizes can be found 
from the law of mass-action, p n = po + n/x+, which re- 
duces to 



p* n = K n (T)pUp+) n e^ +n ^T-^T 



(8) 



where the excess chemical potentials, pf s x = —df el /dp s , 
are 

.PzaLz [ln(f ) + 0(#)] and = 0(jP + ). 

C,(T)/G., as in the 



Rn ex = -i 

h'M'n PT"b 

The equilibrium constant, K(T) 



case of simple electrolyte will remain finite in the limit 
p — > 0. Substituting the expressions for the excess chem- 
ical potential into Eq. || we find that p* w PoP+' ' > 
where g(n) = n — na/T*b + n 2 a/2PT*b. In the limit 
p — > the only possible cluster has the size n* , for which 
the function g(n) attains its minimum. In particular we 
see that for T* < T^ /f — a/b the minimum is attained 
when n* = P(l — T*b/a) counterions are associated with 
a polyion, while for T* > no clusters form and the 
minimum is at n* = 0. In the limit p — > 0, the formation 
of clusters is a continuous transition which occurs when 
the temperature is lowered below T* = = a/b. 
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IV. CONCLUSION 
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As was discussed following Eq. ^ the isomorphism be- 
tween 2d Coulomb law and the interaction potential be- 
tween a polyion and a counterion is valid if b — > 1 (er — > 
—q). In these units Tm — q 2 /ksD, while the Kosterlitz- 
Thouless transition occurs at Tkt — q 2 I^bD, which 
is half the value of the equivalent Manning temperature. 
Furthermore, we would like to stress that while at KT 
transition all the ions associate into dipolar pairs, thus 
leading to a divergent Debye screening length, nothing 
like this happens in the case of polyelectrolytes. Quite on 
the contrary, all the way down to zero temperature there 
remain some free counterions, producing a finite screen- 
ing lenght. Finally, while the KT transition is found at 
non-zero density, the sharpness (discontinuity in slope as 
a function of temperature) of the counterion condensa- 
tion transition will disappear with an increase of density 
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